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r-| . ABSTRACT 

This is a continuation of our study of open clusters based on the 2-Micron All Sky Survey 
photometry. Here we present the results of the mass function analysis for 599 known open clusters 
in the Milky Way. The main goal of this project is a study of the dynamical state of open clusters, 
the mass segregation effect and an estimate of the total mass and the number of cluster members. We 
noticed that the cluster size (core and overall radii) decreases along dynamical evolution of clusters. 
The cluster cores evolve faster than the halo regions and contain proportionally less low-mass stars 
from the beginning of the cluster dynamical evolution. We also noticed, that the star density decreases 
for the larger clusters. Finally, we found an empirical relation describing the exponential decrease of 
the mass function slope with the dynamical evolution of clusters. 

Key words: open clusters and associations: general — infrared: galaxies — astronomical databases: 
2MASS - stars: mass function 

1. Introduction 

The distribution of the stellar masses that stars were formed with, can be de- 
scribed as an empirical function - the Initial Mass Function (IMF, Salpeter 1955, 
Miller & Scalo 1979). According to the basic considerations, the IMF is expected 
to depend on the star forming conditions (Larson 1998). However, the studies pre- 
sented more recently indicate that the mass distribution is relatively invariant from 
one open cluster to another and has a universal character (Phelps & Janes 1993; 
Massey et al. 1995; Durgapal & Pandey 2001; Kroupa 2001, 2002; Pandey et al. 
2005, 2007; Sharma et al. 2008). Similar to our research, Bica & Bonatto (2005) 
and Bonatto & Bica (2005) used the JHK$ photometry data to study the mass func- 
tion (MF) of 5 and 1 1 open clusters, respectively, and pointed out that the MF slope 
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decreases along the dynamical evolution of the investigated clusters. Maciejewski 
& Niedzielski (2007) used homogeneous BV data for 42 open clusters to analyze 
the MF and observed a similar trend. 

This paper is a continuation of our study of known open clusters (Bukowiecki et 
al. 201 1, hereafter Paper I). From our sample of 849 open clusters, we were able to 
analyze the MF for 599 clusters. We used the near-infrared JHK$ photometric data 
to determine the initial mass function (%o)> the dynamical evolution, the mass of 
the cluster and the number of member stars. In Section 2, the method and the data 
analysis is presented. In Section 3, our results are compared to the results found 
in the literature. Section 4 contains discussion of the relations between individual 
parameters. Final conclusions are summarized in Section 5. 

2. Data analysis 

2.1. Data Source and Cluster Selection 

Our research is based on the JHK$ photometric data extracted from the 2MASS0 
Point Source Catalog (Skrutskie et al. 2006). The method is described in Paper I. 
For a sample of 849 open clusters we determined new coordinates of the centers 
and the angular sizes. Moreover, age, reddening, distance, and linear sizes were 
also derived for 754 of them (Table 1 and 2 in Paper I). 

We studied only the MF for a part of the whole sample because of two reasons. 
First, we had to reject clusters without a complete set of parameters. Second, after 
a series of tests, we decided to study clusters for which we observe more than two 
magnitude of the main sequence, between the turnoff point and a level of 15.8 mag 
in the J band - a value of the 99.9% Point Source Catalogue Completeness Limio 
That allowed us to reach a reasonable precision of estimated parameters. After 
applying these criteria, we included 599 open clusters in the further MF analysis. 

2.2. Mass Functions and the Dynamical-Evolution Parameter 

The first step to derive the MF was to build the cluster's luminosity func- 
tion (LF). After a series of tests we used 0.5 mag bins, for the core, the corona 
(halo), and the overall regions separately. Another LF was built for an offset field 
starting at r = r/,„, + 1 arcmin (where r/, m is the limiting size of the open cluster, 
and generally "/ im " is related to the overall cluster). The LF of the offset field was 
subtracted, bin by bin, from the cluster LF, taking the area proportion into account. 
The resulting LF was converted into the MF using a respective isochrone. This way 
we derived the MF approximated by a standard relation of the form 

log(j)(m) = \og{dN /dm) = -(1 +%) \ogm + b Q , (1) 

1 http://www.ipac.caltech.edu/2mass/releases/allsky/ 

2 Following the Level 1 Requirement, according to Explanatory Supple- 
ment to the 2MASS All Sky Data Release and Extended Mission Products 
(http://www.ipac.caltech.edu/2mass/releases/allsky/doc I 
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where m, N , %, and bo are the stellar mass, the number of stars, the mass function 
slope, and a constant, respectively. Using basic astrophysical parameters taken 
from Paper I we estimated the additional cluster parameters: the total mass M tota i , 
total number of stars N tota i , the core mass M core and the number of stars in the core 
Ncore (Fig- !)• Based on our calculations we concluded that 91% of the clusters 
in the sample have a total mass in the range of 100 - 5000 M and the mass of 
the core for 84% clusters is in the range of 50 - 1000 M Q , while 78% of studied 
clusters contain between 100 and 5000 stars. Cores are naturally smaller, 80% of 
the clusters have less than 500 stars within them. These values were calculated by 
extrapolating the MF from the star-mass limit of 0.08M© to the turnoff using the 
methods described in Bica & Bonatto (2005), if % was lower than the universal 
initial mass function, %imf = 1.3 ±0.3 (Kroupa 2001). In the other case the MF 
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Fig. 1. a) Number of clusters in the distribution of the total mass of clusters, b) The same as panel a 
but for core radii, c) Number of clusters in the distribution of the total number of components, d) 
The same as panel c but for core radii. 
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was extrapolated with % = 0.3 from 0.08M Q to the mass of O.5M and for greater 
stellar masses with the derived % to the turnoff. This procedure was applied to 
analysis of the overall MF for the whole sample of 599 open clusters. In addition, 
for 461 clusters, with r core > 0.6 ' , we studied the MF for the core and halo regions 
separately. The error bars for (j)(m) were calculated assuming the Poisson statistics. 
The obtained parameters are listed in Table 1 . 

Section 2 of Paper I describes the procedure of estimating size, age and distance 
in details. 

To describe the dynamic aspect of the studied open clusters, we used the dynamical- 
evolution parameter x , defined as 



where the relaxation time t re i ax was calculated as 



Here t cross =R/o v is the crossing time, N is the number of stars in studied radius 
R and a v is the velocity dispersion (Binney & Tremaine 1987) with a typical value 
of 3 km s _1 (Binney & Merrifield 1998). Our calculations were done for the entire 
region of the open cluster and separately for the core and halo areas. To describe 
differences between the MF slope of the core and the coronal regions we used 
A% = Xhalo — Xcore , which can be treated as the mass segregation indicator. 

2.3. Structural Parameters 

In this section we present the relations of the open cluster's size with the clus- 
ter's mass and the dynamic evolution parameter. From the total sample of studied 
clusters, we analyzed only those for which the relative uncertainty of the total mass 
or the number of potential members was smaller than 100% of the derived value. 
The determinations with greater errors, presented as the open points in all Figures, 
were omitted in further analysis. Analyzing the total mass of 436 clusters (panel a 
of Fig. 2), we obtained a linear relation: 



with the correlation coefficient of 0.70. As one could expect, the core radius is also 
related to the core mass. In this case for 35 1 clusters (panel b of Fig. 2), we fitted 
a linear trend that resulted in an equation: 




cross ■ 



(3) 



\ogR lim = (0.45 ±0.03) log M fofa/ - (0.67 ±0.07) 



(4) 



log/?, 



core 



(0.45 ±0.03) log M t 



core 



(1.13±0.07) 



(5) 



the correlation coefficient was found to be 0.64. These both trends indicate that the 
clusters with small diameters have small mass, and vice versa. A similar conclusion 
was reached by Maciejewski & Niedzielski (2007). 
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Fig. 2. Relations between cluster structural parameters. Full circles represent open clusters used 
for fitting a linear relation and empty circles are determinations with errors greater than 100% of a 
derived value, which have not been taken into account. See Sect. 2.3 for details. 
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We plotted the star density p (stars/parsec 3 ) in panels c and d of Fig. 2 as a 
function of the open cluster's diameter. For the clusters with determined limiting 
(431 objects) and core (350 objects) parameters (the mass and diameter), we found 
linear relations: 

log pa* = (- 1. 80 ± 0.08) log R Um + ( 1.74 ±0.06), (6) 

logp core = (-2. 14 ±0.09) log R core + (1.82 ±0.03) , (7) 

with the correlation coefficients equal to 0.75 and 0.80, respectively. This indicates 
that the star density of open clusters decreases with their diameters. For the smaller 
clusters p/ ;m is much greater than for the larger ones. It is worth noting that for the 
core regions this trend is more significant. This finding suggests that the migration 
of stars in the central regions is stronger. No relation between p and the number of 
cluster members, mass or age was found in the investigated sample. 

The dynamical evolution is described by the x parameter. We used 43 1 (panels 
e and / of Fig. 2) and 350 clusters (panel g of Fig. 2) to investigate how the open 
cluster's core and limiting radii depend on x. We derived the following relations: 

log R Um = (-0. 19 ± 0.02) log xiun + (0.71 ± 0.02) , (8) 

logR core = (-2.20 ±0.02) log %i im + (-0.16 ±0.02), (9) 

\ogR core = (-0.11 ±0.02) log x core + (0.08 ±0.04) , (10) 

resulting in the weak relation with correlation coefficients of 0.53, 0.50 and 0.35, 
respectively. These results demonstrate that open clusters tend to decrease in size 
during their dynamic evolution, both the overall and the core radii. Similar results 
were reported by Maciejewski & Niedzielski (2007) and Nilakshi et al. (2002), 
who noted a decrease in the size for the older and the dynamically more evolved 
systems. 

No relation between the concentration parameter, defined as c = log (Rn m /R C ore) 
(Peterson and King 1975) and x, age or mass segregation was found in the studied 
sample. 

2.4. Mass Function Slopes 

The MF slope is characterized by the % parameter and describes the mass dis- 
tribution in an open cluster. It can be interpreted as a measure of the estimated 
quantities of the low-mass members in comparison to the more massive stars. The 
small x guarantees the deficit of the low-mass stars in the cluster. And, vice 
versa, the higher % illustrate a surplus of low-mass components. We also in- 
spected relations between the MF slope and the dynamical-evolution parameter 
X. Bonatto & Bica (2005) obtained an empirical relation between x and x in the 
form of x(x) = Xo — %i ex P( — W T ) (hereafter the B&B relation; Xo is the ini- 
tial mass function of studied clusters) which suggests the MF slope to decrease 
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Fig. 3. Relations between the MF slopes and other clusters parameters. Full circles represent open 
clusters used for fitting relations %(x) = Xo — Xl ex P (log^Ao) using the solid line and %(x) = %o — 
Xi exp (— Xq/x) with the dotted line (panels a, b and c) or linear relation (panel d). While triangles 
are outlying points and empty circles are measurements with errors greater than 100% of a derived 
value, both (empty circles and triangles) have not been taken into account during fitting functions. 
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exponentially with x. Our results do not confirm this relation. We noticed that 
the MF slope tends to decrease for high values of X. Thus our investigation led 
to another empirical relation describing the exponential decrease more precisely: 
XCO = Xo-Xiexp(logx/x ). 

For Xlim and x, im (panel a of Fig. 3), using the least-square fit to 426 clusters 
and after removing five outlying points, we obtained the following relation: 

XimW = (1.11 ±0.05) - (0.05 ± 0.02) exp (log x ;/m /l. 88 ±0.10) , (11) 

with the correlation coefficients of 0.41 and reduced % 2 = 0.89, while for the B&B 
relation fitted to our data we obtained: 

X Um {%) = (1.09 ±0.05) - (1.31 ±0.24) exp (-52.5 ± 17.6/1,*,) , (12) 

with the correlation coefficients of 0.39 and reduced % 2 = 0.92. 

For Xcore and T core (panel b of Fig. 3), we used 348 clusters and after removing 
two outlying points, we obtained: 

Xcw(x) = (0.86±0.06) - (0.03 ±0.02) exp (log w/0.93 ±0.12) , (13) 

with the correlation coefficient equal to 0.45 and reduced % 2 = 0.88, while for the 
B&B relation, fitted to our data: 

Xcw(x) = (0.82±0.05) - (1.81 ±0.22) exp(-520± 116/w) , (14) 

with the correlation coefficient equal to 0.42 and reduced % 2 = 0.93 . 

For Xhah and X/, m (panel c of Fig. 3), using 328 clusters and after removing 
six outlying points, the least-square fitted relations are: 

XhaioW = (1.19± 0.07) - (0.04 ±0.03) exp (log W0.69± 0.14) , (15) 

with the correlation coefficient equal to 0.36 and reduced % 2 = 0.97 , while for the 
B&B relation: 

Ihaio^) = (1.17 ±0.07) - (1.75 ±0.36) exp (-42.0 ± 14.8/x /im ) , (16) 

with the correlation coefficient equal to 0.34 and reduced % 2 = 1.02. In all these 
cases our new relation fit better than the B&B one. 

It is worth mentioning that the obtained value of %o for both overall and halo 
cluster regions are similar to %imf = 1.3 ±0.3 (Kroupa 2001) while the %q for 
the core is noticeably smaller. The same finding was presented in Maciejewski & 
Niedzielski (2007) and Bonatto & Bica (2005). This suggests that the core regions 
contain proportionally less low-mass stars than the halo areas what may be inter- 
preted as a mass segregation effect. For the core region both relations correlate 
better, indicating that halo regions are treated in some processes that make it diffi- 
cult to study (separating, mixing or errors in estimate of the cluster size). Moreover, 
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the decrease of the MS slope indicates that the evaporation of the low-mass mem- 
bers occurs in the entire cluster. Also the more massive stars might migrate to the 
core from the halo region. 

We noticed weak (the correlation coefficient equal to 0.15) connection between 
the MF slope for overall (panel d of Fig. 3): 

Xum = (0.27 ± 0.08) log(age) + (-1.25 ± 0.62) , (17) 

No relations between the MF slope for core (panel e of Fig. 3) or coronal (panel / 
of Fig. 3) areas and cluster age were detected in the investigated sample. However, 
in the open clusters older than about 500 Myr we noticed that % for all mentioned 
regions becomes either very high or low as compared to the mean value. 

As one can see in panel a of Fig. 4, % core is correlated with the mass segregation 
parameter A%. For 461 open cluster, we obtained the relation: 

%core = (-0.42 ±0.03) A% + (0.74 ±0.06), (18) 

with the correlation coefficient equal to 0.64. As shown in panel b of Fig. 4, Xhalo 
depends on A%, as well. While considering this same number of clusters in the 
previous case, we derived a linear relation: 

Xhalo = (0.58 ± 0.03) A X + (0.74 ± 0.06) , (19) 

with the correlation coefficient of 0.76. 

In panels c, d and e of Fig. 4 we show that the MF slope is not related with the 
location in our Galaxy such as the Galactic Longitude (/), the distance from the 
Galactic plane (Z) or the Galactic center (Rgc)- Panels / and g of Fig. 4 show no 
relations between %lim an d open clusters' structural parameters such as the limiting 
radii (Ru m ) or the concentration parameter (c). Moreover, (panel h of Fig. 4) we 
did not find any relation between the reddening and the MF slope. 

2.5. Mass Segregation 

The dispersion of the MS slope versus the age, similar to that shown in panels 
d,e and / of Fig. 3, can be observed in the comparison presented in panel a of 
Fig. 5. Moreover the mass segregation tends to increase with age, which suggests 
the existence of the initial mass segregation within the protostellar gas cloud, as 
noticed by Maciejewski & Niedzielski (2007) and Sharma et al. (2008). 

The following panel b of Fig. 5 indicates that the mass segregation A% in- 
creases with the dynamic evolution of the core. This in turn suggests that the core 
regions evolve faster than the halo regions so the effect of the evaporation of the 
low-mass members is stronger there. Another process which might be responsible 
for the mentioned increases, is the migration of the more massive components from 
the halo regions to the core. This suggests that the dynamical-evolution parame- 
ter x is a better indicator describing the evolution status of open clusters than the 
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Fig. 4. All panels show relations of the MF slopes with a spatial and structural parameters derived in 
Paper I. See Sect. 2.4 for details. 
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Fig. 5. Evolution of the mass segregation measure in time. 

cluster age. According to equation (2), a small open cluster with small number of 
members can be more evolved than a huge cluster with many components. In panel 
c of Fig. 5 we plotted A% as a function of the dynamical-evolution parameter x . 
No relation can be seen, however only for clusters older than their relaxation time, 
i.e. with log x > we observe a strong mass segregation. A similar results was 
obtained by Maciejewski & Niedzielski (2007). 

In panels a — f of Fig. 6 we present the mass segregation plotted against the 
location in the Galaxy and the structural parameters of open clusters derived in 
Paper I. No statistically relations between these parameters were found. 

3. Comparison with the Published Data 

To test the reliability of our determination, we compared our results with the 
literature data. We collected the literature data from various papers dedicated to 
individual clusters, when available. In panels a and b of Fig. 7 (for 63 and 42 open 
clusters, respectively) presents our MF slope plotted against the published values. 
After rejecting one and six outlying points, respectively, we obtained the linear 
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Fig. 6. All panels show the mass segregation versus location in the Galaxy and structural parameters 
derived in Paper I. See Sect. 2.5 for details. 
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Fig. 7. a) The MF slope for overall derived by us as a function of the literature data, with errors. 
The best-fitting linear relation of the form y = ax is marked with a dashed line. Removed outlying 
points are marked as open circles, b) This same as panel a but for the core MF. c) This same as 
panel a but for the total number of stars, d) This same as panel a but for the total mass. See Sect. 3 
for details. 



relations: 

Xf£ = (0.96±0.Q5)x{L (20) 

XZe = (0-87 ± 0.07) %% re (21) 

with the correlation coefficients of 0.83 and 0.85, respectively. In panel c of Fig. 7 
we present a comparison of the total number of stars in each open clusters. For 44 
objects we received the linear relation: 

NZ a i = (0-94 ±0.09) N% tal (22) 

with the correlation coefficient equal to 0.70. Panel d of Fig. 7 shows a comparison 
of the total mass of each open clusters in the solar masses. For 46 objects we 
obtained the linear relation: 



MZ a i = (0.96 ± 0.09) M'l, 



(23) 
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with the correlation coefficient of 0.66. This illustrates that our analysis of the MF 
is reliable. 

4. Parameter uncertainties 

We discussed the derived parameter uncertainties in section 4.2 of Paper I. Ad- 
ditionally, analyzing the MF, we noticed that the precision of the obtained % param- 
eter depends on the depth of data span between the turnoff and the Completeness 
Limit of the 2MASS photometry (15.8 mag in the J band). We analyzed open 
clusters where this depth was at least 2 mag, but for some of these objects our mea- 
surements exhibit relative uncertainties greater than 100%. We rejected such cases 
from further consideration (open circles in all Figures). An interval of 3 mag range 
or higher would allow to increase the precision of a determinations. However, lim- 
iting over considerations to such clusters would lead to a significant decrease in the 
number of objects. 

While determining the relaxation time we assumed a typical value of the veloc- 
ity dispersion equal to 3 km s _1 (Binney & Merrifield 1998). For some clusters o v 
may have a different value, for example in M67 a v = 0.81 km s _1 according to 
Girard et al. (1989), in Hyades a v = 0.3 km s _1 (Makarov et al. 2000). But even 
the change of velocity dispersion by an order of magnitude, causes only a slight 
change of the log x . 

In panels a and b of Fig. 8 the MF slope and the mass segregation is plotted 
against the distance from the Sun. As one can see, there are no selection effects in 
our sample. Hence our study on the MF seems to be reliable. It is worth mentioning 
that after comparing Xlim and A% with the other parameters derived in Paper I, we 
did not notice any relation effects in the sample neither. 
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Fig. 8. a) Relation between the MF slope and distance from the Sun. b) The same but for the mass 
segregation parameter. See Sect. 4 for details. 
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5. Conclusions 

In this paper, we analyzed the mass function, determined the initial mass func- 
tion (%o)> the mass function slopes (%), the dynamical-evolution parameter x, the 
mass of a cluster and the number of member stars for 599 objects. In the case of 461 
clusters these parameters were derived for cores and the coronal regions separately 

The investigation of our sample leads to the following statistical conclusions: 

o The open clusters with small diameters do not have large masses, and vice 
versa. The star density in open clusters decreases with the cluster's diameter. 

o An average open cluster reach of a medium diameter is less evolved then 
a small cluster with a low number of components. Thus, the dynamical- 
evolution parameter x describes better evolution any status of the open clus- 
ters than just the cluster age. 

o The size of a cluster (both core and limiting) tends to decrease with the 
dynamic-evolution parameter x . 

o The core regions of the open clusters contains proportionally less low-mass 
stars then the halo areas at all cluster ages. This indicates the existence of the 
initial mass segregation. 

o Generally the core regions evolve faster than the halo regions so the effect of 
evaporation of the low-mass members is stronger there. 
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Table 1 

Astrophysical parameters obtained from the mass function analysis. 
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-0.28±6.48 


0.67 ± 3.20 


-2.15 ±7.53 


5 


2 


538 ± 21548 


226 ± 8744 


303 ± 1346 


114 ±498 


King 13 


1.29±0.56 


1.38 ± 1.21 


1.23 ±0.70 


16 


2.17 


13095 ± 5391 


3870 ± 1815 


2675 ±2911 


979 ± 1044 


Berkeley 60 


-0.18±2.28 


2.75 ± 1.03 


-0.78 ± 0.96 





3.93 


702 ± 9027 


449 ± 5770 


27 ±60 


79 ± 178 


Mayer 1 


1.21±0.79 


0.81 ±0.41 


1.10 ±0.51 





6.21 


713 ± 844 


324 ± 383 


279 ± 145 


147 ± 76 


King 1 


2.58±0.91 


1.24 ±0.73 


3.10± 1.10 


13 


1.51 


6620 ± 1649 


2127 ± 526 


1353 ± 263 


461 ± 87 


NGC 103 


1.35±1.36 


0.71 ±0.92 


1.70 ± 1.53 





4.79 


2785 ± 10844 


1164 ±4531 


375 ± 842 


188 ±421 


Stock 20 


0.93±0.55 


0.17 ±0.36 


1.00 ±0.71 





5.06 


867 ± 650 


412 ± 309 


58 ±23 


58 ±23 


NGC 129 


1.02±0.18 


0.43 ± 0.33 


1.02 ±0.17 


1 


5.63 


2964 ± 682 


1408 ± 322 


136 ± 60 


102 ± 45 


Stock 21 


7.82±1.28 


1.60 ± 1.26 


8.66 ± 1.30 


4 


1.69 


1980 ± 972 


601 ± 292 


224 ± 119 


77 ±40 


NGC 133 


0.52±0.43 









11.79 


235 ± 127 


218 ± 117 






NGC 136 


1.18±1.30 


0.48 ±1.82 


1.84 ± 1.12 


4 


2.7 


3526 ±5413 


1365 ±2081 


274 ± 842 


119 ±359 


King 14 


1.23±0.29 


0.92 ± 0.22 


1.15 ±0.36 





8.64 


5901 ± 3117 


2812 ± 1485 


721 ±251 


396 ± 138 


King 15 


-1.42±2.46 


-3.25 ± 2.98 


-1.66 ±2.74 


3 


3.41 


305 ± 2741 


300 ± 2641 


11 ± 130 


28 ± 349 


NGC 146 


1.28±0.27 


0.51 ±0.39 


1.42 ±0.30 





6.25 


5214 ± 1992 


2315 ±885 


382 ± 242 


268 ± 170 


NGC 189 


0.75±0.36 









14.32 


211 ±99 


150 ±71 






Stock 24 


0.82±0.19 


0.13 ±0.43 


0.87 ±0.18 





5.97 


1440 ± 352 


752 ± 184 


54 ±31 


65 ±38 


Dias 1 


0.92±0.50 


1.20 ±0.55 


0.70 ± 0.43 





9.01 


726 ± 658 


404 ± 366 


301 ± 203 


146 ± 99 


King 16 


1.14±0.29 


0.64 ± 0.32 


0.87 ± 0.23 





8.6 


1707 ± 875 


842 ± 432 


266 ± 143 


175 ± 95 


Czernik 2 


0.80±0.33 


0.80 ± 0.53 


0.80 ± 0.27 





15.82 


978 ±416 


685 ± 292 


141 ± 73 


99 ±51 


Berkeley 4 


0.94±0.25 


1.78 ±0.87 


1.33 ±0.39 





21.87 


1415 ±688 


945 ± 459 


289 ± 279 


120 ± 116 


NGC 188 


1.01±0.30 


1.44 ±0.20 


0.74 ± 0.20 


18 


1.37 


12200 ± 597 


3902 ± 190 


3084 ± 100 


992 ± 32 



2 Full data tables available at http://www.astri.uni.torun.pl/- gm/OCS/2mass.html 



